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Abstract 

In the mean field approximation, we study the effects of weak magnetism 
and pseudoscalar interaction in the neutrino energy losses caused by the di- 
rect Urea processes on relativistic nucleons in the degenerate baryon matter. 
Our formula for the neutrino energy losses incorporates the effects of nucleon 
recoil, parity violation, weak magnetism, and pseudoscalar interaction. For 
numerical testing of our formula, we use a self-consistent relativistic model of 
the multicomponent baryon matter. We found that, due to weak magnetism 
effects, relativistic emissivities increase by approximately 40-50%, while the 
pseudoscalar interaction only slightly suppresses the energy losses, approxi- 
mately by 5%. 

PACS number(s): 97.60.Jd , 21.65+f , 95.30.Cq 
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The direct Urea processes, n — > p + l + Ui, p + I — > n + ui, where / is either an electron 
or a muon, are the most powerful neutrino reactions by which neutron stars lose their 
energy. In spite of widely accepted importance of these reactions, the relevant neutrino 
energy losses are not well investigated yet. By modern scenarios, the central density of 
the star can be up to eight times larger than the nuclear saturation density, what implies 
a substantially relativistic motion of nucleons 0. The relevant equation of state of the 
matter is usually calculated in the relativistic approach O, while the energy losses are still 
calculated by the non-relativistic formula obtained by Lattimer et al. [[| more than ten years 
ago. Some aspects of this problem was studied by Leinson and Perez || , who have estimated 
relativistic effects of baryon recoil and parity violation in the direct Urea processes. The 
above relativistic effects increase substantially the neutrino emissivity in comparison with 
the known non-relativistic prediction. 

In the present Letter we derive a relativistic expression for neutrino energy losses caused 
by direct Urea processes on nucleons by taking into account also the effects of weak mag- 
netism and pseudoscalar interaction. It is known that weak magnetism plays an important 
role in core collapse supernovae by increasing mean free paths of antineutrinos ||, |J. Our 
goal is to study the role of weak magnetism in radiation of neutrinos and antineutrinos at 
the long-cooling epoch of neutron star. 

We employ the Walecka-type relativistic model of baryon matter J7|], where the baryons 
interact via exchange of a, lu, and p mesons, and perform the calculation of the neutrino 
energy losses in the mean field approximation. This approximation is widely used in the 
theory of relativistic nuclear matter, and allows to calculate in a self-consistent way the 
composition of the matter together with energies, and effective masses of the baryons. The 
Lagrangian density, which includes the interaction of a nucleon field \I/ with a scalar field a, 
a vector field w„ and an isovector field b„ of p-meson is of the formBi 



£ = * 



- ( id" - 9^ - ^gpnW ■ t ) - ( \ [ u - < hD <7 ) 



+ l - (d.ad^a - m 2 y) - U (a) + l(i lfi d» - m,) I, (1) 

Here ^ are the Dirac spinor fields for nucleons, b M is the isovector field of p-meson. We 
denote as r the isospin operator, which acts on the nucleons of the bare mass M. The 
leptons are represented only by electrons and muons, I = e~,[i~, which are included in the 
model as noninteracting particles. The field strength tensors for the uj and p mesons are 
Fpu = dfj,uj v — dyUJ^ and = d fJi b u — <9„b M , respectively. The potential U (a) represents the 
self-interaction of the scalar field and is taken to be of the form 



In principle, the pion fields should be also included in the model. However, the expectation 
value of the pion field equals zero, giving no contribution to the mean fields. Therefore, only 
non-redundant terms are exhibited in the Lagrangian density. 

2 In what follows we use the system of units h = c = 1 and the Boltzmann constant kg = 1. 
Summation over repeated Greek indexes is assumed 
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U (a) = -bM {g aN af + -c (g uN a) 4 . (2) 

In what follows we consider the mean field approximation widely used in the theory of 
relativistic nuclear matter. In this approximation, the meson fields are replaced with their 
expectation values 

<7 - (<7> = <7 , /^KIewo^, b"-(b") = (0,0,A))<W (3) 

In this case only the nucleon fields must be quantized. This procedure yields the following 
linear Dirac equation for the nucleon 

Udp-f - g u -fu - -9 P l°p r 3 — (M — g a a )^j * (x) = 0, (4) 

Here and below we denote as T3, and r± = {r% ±2x2) /2 the components of isospin operator, 
which act on the isobaric doublet \l/ (x) of nucleon field. 

The stationary and uniform condensate fields equally shift the effective masses 

M* = M-g a a , (5) 

but lead to different potential energies of the proton and neutron 

U n = gjjj - ^g p p , U p = g^uo + ^g p p , (6) 

thus creating the energy gap U n — U p = —g p Po between possible energies of protons and 
neutrons. 

The exact solutions of Eq. (f|) can be found separately for protons and for neutrons. In 
our case of a stationary and uniform system, solutions are the spinor plane waves 

V> n (z) = N n u Q (P) exp (-iE n t + ipr) , (7) 



% (x) = N p u p (P') exp {-iE p t + ip'r) , (8) 

where the neutron and proton energies are given by E n (p) = \/p' 2 + M* 2 + U n , and E p (p') = 
\/p 12 + M* 2 + U p . The free-like spinors u n (P) and u p (P') are constructed from the kinetic 
momenta of the neutron and the proton 

= (E Q - U n , p) = (v/p2 + M* 2 , p) , (9) 



P"> = (P p - U p , p') = (^W, p') . (10) 

In what follows we denote by e = a/p 2 + M* 2 , e' = y/p' 2 + M* 2 the kinetic energy of the 
neutron and the proton respectively. So the normalization factors are 

Nn = -|U N p = -=, (11) 
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and the single-particle energies can be written as E n (p) = e + U n , and E p (p') = e' + U p . 

We consider massless neutrinos of energy and momentum ki = (o;i,ki) with ui = |ki| . 
The energy-momentum of the final lepton I = e~, fi~ of mass mi is denoted as k 2 = (u 2 , k 2 ) 



with to 2 = y k 2 + mf. In the lowest order in the Fermi weak coupling constant Gp, the 
matrix element of the neutron beta decay is of the form 

(f\ (S - 1) \i) = -i^N n N p u t (k 2 ) 7/4 (1 + 75) v {-k x ) p (P'\ J» (0) |P) n x 

x (2vr) 4 5 {E B -E v -u x - u 2 ) 8 (p - p' - k x - k 2 ) , (12) 

were C = cos 9c = 0. 973 is the Cabibbo factor. The effective charged weak current in the 
medium consists of the polar vector and the axial vector, J 71 (x) = V^ 1 (x) + (x) . 

Our goal now is to derive the nucleon matrix element of the charged weak current in the 
medium. Consider first the polar-vector contribution. The Lagrangian density (|l|) ensures 
a conserved isovector current [§] 

X m = l^rif; + h u x B"", dfjT^ = 0. (13) 

Besides the directly nucleon contribution this current includes also the contribution of the 
isovector field b M , which obeys the field equations 

d v W* + m 2 y = ^g p $'fTil>, d v h v = 0. (14) 
By the use of Eq. flT4]) the condition c^T^ = may be transformed as 

id p (^7^) = g p b p x $7^. (15) 
In the mean field approximation this gives 

id P (f'fT+i/^ = -gpPo^T+ii. (16) 
By introducing the covariant derivative 

(^-ig pPo ,v\, (17) 



we can recast the Eq. (16) to the following form 

Dp (^r+V) = 0. (18) 
At the level of matrix elements this can be written as 

u p (P') q^fu* (P) = 0, (19) 
where q p is the kinetic momentum transfer 

<f = (E n -E p + g pPo , p - p ; ) = (e - e', p - p') (20) 
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Thus the matrix element of the transition current is orthogonal to the kinetic momentum 
transfer, but not to the total momentum transfer from the nucleoni. Note that this effect 
originates not from a special form (Jf^) of the conserved isovector current in the medium 
but is caused by the energy gap between the proton and neutron spectrums. This follows 
directly from the Dirac equation ([J), which ensures the Eq. ([19]) with 

t = ~ P'* 1 = (En -E p -U n + U p ,p- p') , (21) 

which coincides with Eq. ( pOD because U n — U p = —g p po- 

By the use of the isovector current (0) one can construct the conserved electromagnetic 
current in the medium 

JL = + n + ^7 d » (*Aa^*) , <Ve M m = 0. (22) 

The last term in Eq. ( |2~2|) is the Pauli contribution, where 2a fJlU = — and 

A = A p ^(l + r3) + A n i(l-r 3 ). (23) 

In the mean field approximation, we replace the magnetic formfactors of the nucleon with 
anomalous magnetic moments of the proton and the neutron, A p = 1.7928 and A n = —1.9132. 

By the conserved- vector-current theory (CVC), the nucleon matrix element of the charged 
vector weak current is given by 

p (P'\ V» \P) n = p (P'\ Jt? m \P) p - n (P'\ J» m \P) Q . (24) 

Thus, in the mean field approximation, we obtain 



,<P'| V» (0) \P) n = u p (P') 



7 + 2M Qv 



u n (P) . (25) 



The second term in Eq. (p5[), describes the weak magnetism effects. By the use of Eq. (|J) 
we find 

P (P'| (0) \P) n = 0. (26) 

Consider now the axial- vector charged current. This current is responsible for both the 
np transitions and the pion decay. In the limit of chiral symmetry, m n — > 0, the axial-vector 
current must be conserved. In the medium with p meson condensate, this implies 

lim D^A" (x) = 0, (27) 

m-„-— >0 

where the covariant derivative is defined by Eq. (]TTf) - At the finite mass of a pion, 
m n , the axial- vector charged current is connected to the field 7r_ (x) = (iti + in 2 ) /v2 of n~ 



3 Some consequences of this fact for the weak response functions of the medium are discussed in 
1 
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meson. For a free space, this relation is known as the hypothesis of partial conservation of 
the axial current (PCAC). In the medium the PCAC takes the form 

(x) = mlUir- (x) , (28) 

where = 139 MeV is the mass of 7r-meson, and f n is the pion decay constant. 

With allowing for interactions of the pions with nucleons and p mesons the Lagrangian 
density for the pion field is of the form H 
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(d^TT - g p b^ x ?r) ■ (d^ir - g^ x vr) - m^ir • tt + ig^^r ■ Trip. (29) 



In the mean field approximation this results in the following equation for the field of tt 
meson 

[JidP + g p p ) 2 - (iV) 2 - mlj vr_ (x) = -V2ig v ijj p <y 5 iJ) lx , (30) 

where g n is the pion-nucleon coupling constant. 

For the nucleon transition of our interest, the Eq. (|28|) gives 



p (P'\ q„A» (0) \P) n = imlU p (P'\ tt (0) |P) n (31) 
Here the right-hand side can be calculated by the use of Eq. fl30|). We obtain 

p (P'| (0) \P) n = ~ ^ l !f u p (P') 75Mn (P) . (32) 

This equation allows to derive the nucleon matrix element of the axial- vector charged current. 
Really, to construct the axial- vector matrix element of the charged current, caused by the 
nucleon transition, we have only two independent pseudovectors, consistent with invariance 
of strong interactions under T 2 isospin transformation, namely: u p (P') 7 M 7 5 u n (P), and 
u p (P') g Ai 75M n (P). This means that the matrix element of the axial- vector charged current 
is of following general form 

P (P'\ A» (0) |P) n = C A u p (P') (7^75 + F q q^ 5 ) u n (P) . (33) 

Here, in the mean field approximation, we set C A — 1.26, while F q is the form-factor to be 
chosen to satisfy the Eq. (0), which now reads 

C A (-2M* + F q q 2 ) u p (P') 75 «n (P) = - ^^t^ P ( P ') 75«n (P) . (34) 

m n q 

Thus 

C A (2M* - F q q 2 ) = ^Mll^L. (35) 

In the mean field approximation, we assume that the coupling constants are independent of 
the momentum transfer. By setting q 2 = in Eq. fl35|) we obtain the Goldberger - Treiman 
relation f n g n = \/2M*C A . By inserting this in ( |35l) we find 
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2M* 



(ml - q 2 ) 



Thus, with taking into account Eqs. 
neutron beta decay is found to be 



M 



fi 



GpC 

-i — — ui (k 2 ) 7 M (1 + 7s) v {-h) x 



V2 



(36) 

3D, and (|36"D , the total matrix element of the 

(37) 



CvY + ^G M a^q v + C A (7^75 + F q q^ 5 ) 



xm p (P' 

where, in the mean field approximation, we assume 



u a (P) 



c 



M 



Xp — X r 



3.7, C A = 1-26. 



(38) 



Note that the matrix element obtained is of the same form as that for the neutron decay 
in a free space, but with the total momentum transfer replaced with the kinetic momentum 
transfer. Due to the difference in the neutron and proton potential energy, the kinetic 
momentum transfer 



p-p> 



(£ - £ , p - p 



(39) 



to be used in the matrix element (|37|) differs from the total momentum of the final lepton 
pair 



K = k x + k 2 = (e - e' + U n - U p , p - p') 



(40) 



This ensures K 2 > 0, while q 2 = (e — e') 2 — (p — p') 2 < . 

The square of the matrix element of the reaction summed over spins of initial and final 
particles is found to bei: 



\M 



fi\ 



32G|C 2 



(C\ - C 2 V ) M* 2 (hk 2 ) + (C A - C v f (hP 2 ) (hPt) 



+ (C A + C v ) 2 (k 1 P 1 )(k 2 P 2 ) 

+Cm 1T [2Ca ^ klPl) ^ hP2) ~ ^ klP2) ^ Pl)) 

+C V {{k x k 2 ) (PxP 2 - M* 2 ) - (k x Pi - hP 2 ) (k 2 P, - k 2 P 2 



AM 2 



M* (hP 2 )(3 (k 2 P 2 ) - (hPi)) 



+M* 2 (hPi) (3 (k 2 Pi) - (k 2 P 2 )) + (hk 2 ) (PiP 2 - M* 2 

- (hPx + hP 2 ) (k 2 Px + k 2 P 2 ) (P X P 2 )] 

+C 2 A F q (2M* + F q (M* 2 - (PP 2 ))) [(hk 2 ) (M* 2 - (P X P 2 )) 

- (hPi - hP 2 ) [k 2 P x - k 2 P 2 )\] 



(41) 



4 Here we rectify an error made in the journal version of the paper. The correct expression can 

by simple replacement Cm 



be obtained from that published in [10|, |T 
grateful to M. Prakash and S. Ratkovic who have pointed out this error 



Cm/2. The author is 
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with Pi = (e, p) and P 2 = (e', p'). 

We consider the total energy which is emitted into neutrino and antineutrino per unit 
volume and time. Within beta equilibrium, the inverse reaction p + l — > n + ui corresponding 
to a capture of the lepton /, gives the same emissivity as the beta decay, but in neutrinos. 
Thus, the total energy loss Q for the Urea processes is twice more than that caused by the 
beta decay. Taking this into account by Fermi's "golden" rule we have 

f d 3 k 2 d 3 k l d 3 pd 3 p' 2 , /, £ W1 ,x 

x (2tt) 4 S (E„ (p) - E p (p') -u l - wj) S (p - p' - Id - k 2 ) . (42) 

Antineutrinos are assumed to be freely escaping. The distribution function of initial neutrons 
as well as blocking of final states of the proton and the lepton I are taken into account by the 
Pauli blocking-factor f n (1 — f p ) (1 — fi). The Fermi-Dirac distribution function of leptons 
is given by 

fi M = - ( 1 WT ■ n . (43) 

exp (u> 2 - lii) /T+l 

while the individual Fermi distributions of nucleons are of the form 

/n {£) = ex P ((e + U n -fi n )/T) + V (44) 



Ip[£) exp((e' + U p -^)/T) + V 

By neglecting the chemical potential of escaping neutrinos, we can write the condition of 
chemical equilibrium as p,i = p> n — fi p . Then by the use of the energy conservation equation, 
£ + U n = e' + Up + 0J2 + uji, and taking the total energy of the final lepton and antineutrino 
as 0J2 + oj\ = Hi + uj' we can recast the blocking-factor as 

/ n (a) (1 - / p (cO) (1 - fi W) 
= fn (e) (1 -fn(e- u')) (1 - /, {m + u f - , (46) 

where uj' ~ T. 

Furthermore, since the antineutrino energy is U\ ~ T, and the antineutrino momentum 
|ki| ~ T is much smaller than the momenta of other particles, we can neglect the neutrino 
contributions in the energy-momentum conserving delta-functions 

S(e + U n -e' -Up-ux- uj 2 ) 5 (p - p' - k x - k 2 ) 
~ 5 (e + U n - e' - Up - uj 2 ) 6{p-p'- k 2 ) (47) 

and perform integral over d 3 p' to obtain p' = p — k 2 in the next integrals. 

The energy exchange in the matter goes naturally on the temperature scale ~ T, which is 
small compared to typical kinetic energies of degenerate particles. Therefore, in all smooth 
functions under the integral (|42|) , the momenta of in- medium fermions can be fixed at their 
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values at Fermi surfaces, which we denote as p n , p p for the nucleons and pi for leptons 
respectively. 

The energy of the final lepton is close to its Fermi energy /// = // n — /i p , and the chemical 
potentials of nucleons can be approximated by their individual Fermi energies /i n = e n + U n , 
Hp — £p + Up. This allows us to transform the energy-conserving 5-function as 



<^n - \/Pn + P\ ~ ZPnPl COS 9 X + M* 2 + U n - U p - /I X ) 



^h'-^M-^)' (48) 

where 9 X is the angle between the momentum p n of the initial neutron and the momentum 
p x of the final lepton. Notice, when the baryon and lepton momenta are at their individual 
Fermi surfaces, the 5- function fl48[ ) does not vanish only if p p + pi > p n - 

Further we use the particular frame with Z-axis directed along the neutron momentum 
p n . Then 

F a = (0, 0, p n , e n ) 

k\ = lo\ (sin^, 0, cos 6^, 1) 

h = (pi sin 6 1 cos <p h pi sin 8 X sin ip x , p x cos 8 U , fi x ) (49) 

The energy-momentum of the final proton is defined by conservation laws: 

P2 = {~Pi sin 81 cos (pi, -pi sin 9 X sin (p t , p D - p x cos 9 V , e p ) (50) 

Insertion of ( pE9|) and (|5T| ) in the square of the matrix element fl4T|) yields a rather cumbersome 
expression, which, however, is readily integrable over solid angles of the particles. 

Since we focus on the actually important case of degenerate nucleons and leptons, we 
may consider the neutrino energy losses to the lowest accuracy in T/fi x . Then the remaining 
integration reduces to the factor 

duj x uj\dJ de f n (e) (1 -f n {e- oj')) (1 - /, + J - cut)) 
00 , , J r 00 , cj? 457 



dJ / duj x — - = ^^7T 6 T 6 . (51) 

exp u'/T - 1 Jo 1 + exp (wi - w') /T 5040 v ; 



Finally the neutrino emissivity is found to be of the formi: 

Q = if^^ (Pi + Vp - Pn) { {C\ - C v ) M* 2 fi x 



2 2 
-Pi 



+\ {pv + c a) 4e n £ P /i« - (e n - £ P ) ((e n + e p ] 
A/[* 

+CvCm 2M I 2 ^ n ~ £p ^ ^ ~ ( 3 ^ n ~ ~ P ^ ^ 



5 See footnote to Eq. @ 
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+C A ( C v + -tjCm J (en + £ P ) (pf - (e n - e 



? 

V) 



M 

+ (p 2 -(e n - £p ) 2 ) (24 + 24-^) /i/ 

- (pf - (e n - £ p ) 2 ) (e„ + £ p ) 2 (2e: n - 2e p - m) 

-CiM* 2 $ (l + m 2 $) [/u ((e n - e p f + pj) - 2 (e n - e p ) p 2 ] } (52) 
with (2) = 1 if x > and zero otherwise. In the above, the last term, with 

$ = 2x 2 \ 3' ( 53 ) 

m l + Pi ~ On - £p) 

represents the contribution of the pseudoscalar interaction. The "triangle" condition p p + 
pi > p n , required by the step-function, is necessary for conservation of the total momentum 
in the reaction and exhibits the threshold dependence on the proton concentration. 

Neutrino energy losses caused by the direct Urea on nucleons depend essentially on the 
composition of beta-stable nuclear matter. Therefore, in order to estimate the relativistic 
effects, we consider the model of nuclear matter, which besides nucleons includes E and A hy- 
perons |J. The parameters of the model are chosen as suggested in Ref. to reproduce the 
nuclear matter equilibrium density, the binding energy per nucleon, the symmetry energy, 
the compression modulus, and the nucleon effective mass at saturation density n = 0.16 
fm~ 3 . The composition of neutrino- free matter in beta equilibrium among nucleons, hy- 
perons, electrons and muons is shown on the left panel of Fig. 1 versus the baryon number 
density n b , in units of n . 

On the right panel of Fig. 1, by a solid curve we show the relativistic neutrino emissivity 



of reactions n — > p + e~ + u e , p + e~ — >• n + u e , as given by Eq. (|52|) . Appearance of 
hyperons in the system suppresses the nucleon fractions and lepton abundance. Therefore 
at densities, where the number of hyperons is comparable with the number of protons, the 
relativistic emissivity reaches the maximum and then has a tendency to decrease. To inspect 
the contribution of weak magnetism and the pseudoscalar interaction we demonstrate two 
additional graphs. The long-dashed curve demonstrates the energy losses obtained from Eq. 
( |52D by formal setting $ = 0. This eliminates the pseudoscalar contribution. The dot-dashed 
curve is obtained by formal replacing $ = and Cm — 0, which eliminates both the weak 
magnetism and pseudoscalar contributions. A comparison of these curves demonstrates the 
weak magnetism effects. The contribution of the pseudoscalar interaction can be observed 
by comparing the total neutrino energy losses (solid curve) with the long-dashed curve, 
which is calculated without this contribution. We see that, due to weak magnetism effects, 
relativistic emissivities increase by approximately 40-50%, while the pseudoscalar interaction 
only slightly suppresses the energy losses, approximately by 5%. 
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Fig. 1. The left panel shows individual concentrations for matter in beta equilibrium 
among nucleons, hyperons, electrons and muons as a function of the density ratio nh/n Q . 
The right panel represents the neutrino emissivity of the direct Urea processes among nu- 
cleons and electrons for the matter composition represented on the left panel. The curves 
begin at the threshold density. The solid curve represents the total relativistic emissivity, 
as given by Eq. ([52]). The dot-dashed curve shows the relativistic emissivity without con- 
tributions of weak magnetism and pseudoscalar interaction, and the long-dashed curve is 
the emissivity without the pseudoscalar contribution. All the emissivities are given in units 
10 27 Tg ergcm _3 s _1 , where the temperature T 9 = T/10 9 K. 



12 



